Suppression of d-wave superconductivity in the checkerboard Hubbard model 
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Using a dynamical cluster quantum Monte Carlo approximation we investigate the d-wave su- 
perconducting transition temperature T c in the doped 2D repulsive Hubbard model with a weak 
inhomogeneity. The inhomogeneity is introduced in the hoppings t' and t in the form of a checker- 
board pattern where t is the hopping within a 2 x 2 plaquette and t' is the hopping between the 
plaquettes. We find inhomogeneity suppresses T c . The characteristic spin excitation energy and the 
strength of d-wave pairing interaction decrease with decreasing T c suggesting a strong correlation 
between these quantities. 



Introduction. The role of inhomogeneity in High Tem- 
perature Superconductors (HTS) is still an unsettled is- 
sue. Neutron scattering experiments reveal the pres- 
ence of ID charge and spin microscopic inhomogencities 
(stripes) in under-doped HTS 0,0,13,11 above the super- 
conducting transition temperature T c . Scanning tunnel- 
ing microscopic imaging provide evidence for the presence 
of two dimensional real space modulations or local den- 
sity of states modulations called checkerboard patterns 
in Bi-2212 0,0,0,(1 and Na-CCOC @, [Io| supercon- 
ducting samples. Inhomogencities observed in cuprates 
led to theoretical scenarios for an inhomogeneity based 
pairing mechanism in HTS [ll|, 12 1. Other theoretical 
studies!!; 
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17] argue that inhomogencities can 
enhance pairing and T c . 

The checkerboard model[l7| we investigate, is sketched 
in the upper panel of Fig. [TJ It is homogeneous when 
t' = t and inhomogeneous when t' / t. Here t(t') is 
the hopping amplitude within (between) the cluster(s). 
The low-temperature phase diagram of the checkerboard 
Hubbard model (CBHM) in the strongly inhomogeneous 
limit (t 1 <C t), i.e with weakly coupled clusters, was stud- 
ied in Refs. [l^, 0. For small t' <i and for values of 
on-site interaction U ~ 4.58t, the authors find a d-wave 
superconducting state in addition to a variety of other 
phases. Their results led them to expect a maximum of 
T c at an optimal t' < t. Other authors [13, (2(| find d- 
wave superconducting state with a significant T c for the 
homogeneous case (t' = t) in the 2D Hubbard model. 
Therefore, it is interesting to see if the inhomogeneous 
CBHM has a higher T c compared to the homogeneous 
model. This is especially important in the weak inhomo- 
geneity regime not addressed before, since experimentally 
observed inhomogencities in the cuprates are weak. 

In this paper we investigate <i-wave superconductiv- 
ity in the CBHM near the homogeneous limit, i.e., for 
a weakly inhomogeneous system. We use the dynamical 
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Figure 1: (color on-line) 2x2 Plaquette model. Here t is 
intra-cluster hopping and t' is inter-cluster hopping. 



cluster approximation (DCA) [24[ to calculate T c as a 
function of t' . DCA is a momentuni(fc) space formula- 
tion and suitable to study the problem since our interest 
is around the homogeneous limit. We find that a weak 
inhomogeneity suppresses T c . At fixed values of U/W 
and U/t with doping S appropriate for HTS, T c decreases 
for t' =/= t. Therefore, we find a maximum T c in the ho- 
mogeneous system. Furthermore calculations with fixed 
U/t for t' < t, show that T c decreases proportional to the 
characteristic spin excitation energy and the strength of 
the d-wave pairing interaction, suggesting a strong cor- 
relation between them. 

Formalism. The DCA maps the lattice problem onto 
a periodic cluster embedded in a self consistent host. 
Short range correlations up to the linear cluster size L c 
are treated explicitly while longer range correlations are 
treated in a dynamical mean-field manner. We use a 
Quantum Monte Carlo (QMC) method to solve the clus- 
ter problem. We use the Maximum Entropy method to 
analytically continue imaginary time QMC data to real 
frequencies [26j |. 

The Hamiltonian of our model is 



(ij)a i 

Where c\ a {cia) creates (destroys) an electron at site i 
with spin a and ni a = c\ a Ci a . Here (ij) denotes the 
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nearest neighbor sites i and j. As shown in Fig. [TJ t%j = 
t(tij = t'), when i and j belong to the same (neighboring) 
plaquette(s). U is the local Coulomb repulsion and it is 
site independent. 

Here we present results obtained for a N c = 4 sites 
cluster (2x2 sites). N c = 4 is the smallest cluster which 
contains one d-wave plaquette and allows for a d-wave 
superconducting state [19||. It has been intensively stud- 
ied in the last few years |25|]. It has the advantage that 
the QMC sign problem is mild thus allowing the investi- 
gation of physics at low temperatures. However, because 
it contains only one d-wave plaquette it does not cap- 
ture d-wave phase fluctuations and hence overestimates 
T c [l^]. Calculations for larger clusters that capture these 
fluctuations will be an interesting future work. Details 
of the full DCA formalism applied to the CBHM will be 
published elsewhere 2l| . 

T c and the Neel temperature (Tjv) are the tempera- 
tures where the superconducting and antiferromagnetic 
susceptibilities diverge. For small clusters it may be 
approximately described as the temperature at which 
the superconducting correlation length exceeds the lin- 
ear cluster size. In a mean field calculation like DCA the 
transition temperature is overestimated and decreases 
with increasing cluster size [25I ]. As we are treating 
the model in an elaborate mean-field scheme, the phase 
transition signalled by the divergence will quite likely 
not be present in the real 2D system, but may sur- 
vive as Kosterlitz-Thouless transition, at least for d-wave 
superconductivity ■ 

T c can be obtained from the temperature dependence 
of the pairing matrix M(K, K') = T PP {K, K')x°(K', K). 
Here, K = (K,iu n ), T PP (K, K') (pairing interaction) 
is the irreducible particle-particle vertex and the bare 
bubble, x°(K',K), is the coarse-grained product of two 
fully dressed single particle Green's functions G(k, k') 
and G(—k, —k'). ui n = (2n + 1)ttT is the Matsubara 
frequency at temperature T. We determine T c from the 
temperature dependence of the leading eigenvalue A(T) 
of M [H, ■ We find T c from A(T C ) = 1 and the sym- 
metry of the superconducting state is determined by the 
symmetry of the corresponding eigenvector, <fr(K), of M. 
For the investigated range of t' , we find that the leading 
eigenvector &d(K) has the d-wave symmetry. 

We define the strength of the pairing interaction Vd, 
by projecting out the pairing vertex r pp , to the d-wave 
sub-pace as given in Eqn. [2] 



Vd(T) = ^" E 3(K)r^(K,7rT;K',7rT) 5 (K') 



(2) 



KK' 



Here g(K) = cos(K T ) — cos(Kj,) is the d-wave form fac- 
tor. Correspondingly, we define the d-wave projected 
bare bubble 

PMT) = ^ E 2(K)x° p (K, K', ttTMK') (3) 



KK' 



Vd, Pdo were initially introduced in Ref. [29(. They sat- 
isfy the approximation V d (T)P d0 (T) w \ d {T). 

Results. Since the band width W = 4(t + t') of the in- 
homogeneous system changes with t' , we could use either 
W or the maximum of t and t' as a unit of energy. Since 
both the Hamiltonian and W are symmetric under the 
interchange of t and t', in both cases the results should 
reflect this symmetry. The first choice compensates for 
changes in T c due to changes of the kinetic energy. The 
second choice reflects the hopping suppression due to in- 
homogeneity and will produce perfectly symmetric plots 
with abscissa t'/t around one. 

We find that both of these choices lead to the same 
main conclusion that the inhomogeneity suppresses T c . 
In Fig. [2] (a) we show results when W is taken as the 
energy unit, for U/W = 1 and at 10% doping. T c is 
maximum in the homogeneous limit t'/t = 1. In the 
subsequent figures throughout the paper we will present 
results for the case where the energy unit is given by t. 
As can be seen from Fig. (2jb), where U = 6t at 8% dop- 
ing case is shown, T c monotonically decreases with de- 
creasing t' . The suppression of T c due to inhomogeneity 
occurs on the doping range appropriate for HTS. This 
can be seen in Fig. HJc), where we show T c vs. 5 (2?| 
for U = 8t. Note that the inhomogeneous system with 
t 1 = 0.88i has a lower T c than the the homogeneous sys- 
tem on the doping range 2% — 20%. 

We find that t' < t enhances the density of states N(oj) 
at low excitation energy (— 0.5f < u> < 0.5t). This can be 
seen in Fig.^d) where N(w) for t' = 0.88i and t' = t arc 
compared at 8% doping for U = 8t. The enhancement 
of N(uj) is a band-structure effect, since it can also be 
noticed for the non-interacting problem, i.e., when U = 0. 
The low energy part of N(ui) has an important role in 
determining the value of T c . 

We find two competing effects of t' on the pairing ma- 
trix as captured by Vd and Pdo- Fig. [3ja) shows the t' 
dependence of Vd, Pdo and T c for U — 6t at 8% dop- 
ing. Vd(t') and Pdo(t') are shown at a temperature equal 
to the homogeneous T c ; ie. T c -% om , and arc normal- 
ized to their values at t' = t for easy comparison. Vd 
decreases with decreasing t', indicating a decreasing d- 
wave contribution to T pp as t' decreases. The effect of 
if on x° is captured by Pdo- Pdo(t') increases with de- 
creasing t'. Since Pdo is related to N(lu) through G, this 
behavior is related to the enhanced low energy part of 
N(u>) discussed previously. The effect of t' on Vd and 
Pdo provide a qualitative description of the behavior of 
T c since V d {T)P d0 (T) w X d {T). Note that the relative 
decrease in Vd is always larger than the relative increase 
in Pdo- Therefore, the net effect of t' is to reduce T c . 
The behavior of Vd(t') and Pdo(t') at other values of U 
and d (not shown) is similar, leading to a suppression of 
superconductivity. 

Using the DCA for the homogeneous 2D Hubbard 
model, it was previously shown that the dominant part 
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Figure 2: (color on-line) (&)T c {t')/W for U/W = 1 and 10% 
doping. T c is suppressed for t' 7^ t. (b) T c (t')/t for U — 6t 
and 8% doping. T c monotonically decreases with decreasing 
t' . (c) T c /t for t' — t and t' = 0.88t as a function of doping 
S when U — 8t. (d) Low energy density of states N(u)t for 
t' = t and t' = 0.88t at T = O.llt for 8% doping and U = 8t. 
Note the increase in low energy spectral weight for t' = 0.88t. 



of the pairing interaction comes from the 5=1 particle- 
hole magnetic channel [23[. The energy scale of spin ex- 
citations is characterized by the super exchange interac- 
tion, J = At 2 /U and changing the hopping matrix ele- 
ments of the homogeneous problem will also influence J. 
However, the relationship between T c and the character- 
istic energy scale of spin excitations, in the weakly inho- 
mogencous system is not known. Therefore it is worth 
studying the effect of t' on the Neel temperature and 
on the spin excitations as characterized by the magnetic 
structure factor S(Q,lj) 

We find that at small doping, Tjy decreases with de- 
creasing t' as shown for U = St at 2% hole doping in 
Fig. [3](b) . As discussed before, for a finite size cluster 
Tjv occurs when the antiferromagnetic correlation length 
(af exceeds the linear cluster size. Therefore decreasing 
Qaf indicates the suppression of antiferromagnetism in 
the weakly inhomogencous system for small 8 compared 
to a homogeneous system. 

The magnetic structure factor S[Q = (0, 7r), uj] for var- 
ious t' values is shown in Fig. 0Ja) for U = 6t at 8% 
doping. The location of the peak of 5[Q = (0, n),u>] can 
be considered a measure of the characteristic spin excita- 
tion energy and we define its position to be 2J e //. This 
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Figure 3: (color on-line) (a) Normalized values of T c , Va and 
Pdo as a function of t' at 8% doping, U = 6t and N c — 4. Va 
and Pdo are described in the text, (b) Neel temperature Tjv/i 
for different values of t' at 2% doping for U — 8t. 
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Figure 4: (color on-line) (a) Magnetic structure factor S(Q, uj) 
for U = 6t, 6 = 0.08 and temperature T = 0.087£ at Q = (0, tt) 
for different values of t' . The location of the peak is a measure 
of the characteristic spin excitation energy, here defined as 
Jeff - Inset of (a) is a blow up to show the area around peak 
positions. T c (t') and J e //(t',T = 0.087£) normalized to their 
values at t' = t, as a function of i! for (b) U = 6t, S = 0.08, 
(c) U = 8t,S = 0.08 and (d) U = 8t,S = 0.15. Note that both 
normalized T c (t') and Jeffit') monotonically increase with t . 



is done in analogy with linear spin wave theory [3Ct |32| 
where the peak of the structure factor at Q = (0, 71") is 
at energy 2 J, J being the super exchange nearest neigh- 
bor spin interaction. Note that J e ft(t') decreases with 
decreasing t\ showing that a weak inhomogeneity soften 
the high energy spin excitations. This can be seen from 
the plot of 5[Q = (0, 7t),cj] in Fig. 0|a) and plots of 
Jeffit 1 ) in Figs. Hb), (c) and (d). 

We find that T c cx J e ff, as suggested from Figs. HJb), 
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(c) and (d). By comparing 8% doping Figs.[4jb) and (c) 
with 15% doping Fig. BJd), one can see that this propor- 
tionality holds better for small doping. The propotion- 
ality between T c and the super exchange interaction J 
are in accordance with the experimental observations re- 
ported in Ref. [33| and the computational investigation 
on the homogeneous system in Rcf . [2(| . 

Discussion. The maximum T c occurs in the homoge- 
neous system and it is suppressed by weak inhomogenc- 
ity. The suppression of T c prevails over the investigated 
range of doping 8 and U, which are appropriate for a 
description of HTS. A recent experiment on Na doped 
Ca2Cu02Cl2 (Na-CCOC) also indicates suppression of 
T a due to inhomogeneity over a wide range of doping 
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According to our analysis, the suppression of T c due 
to weak inhomogeneity is accompanied by a suppression 
of Vd and J e / / . This is an important observation since it 
suggests a strong correlation between the characteristic 
spin excitation energy J e // and the pairing interaction. 
It is consistent with the previous finding that the pairing 
interaction in the 2D Hubbard model is dominated by the 
5=1 magnetic channel |23j and other results reported 
by Maier et al. [1^, . Therefore inhomogeneities that 
reduce magnetic contributions to the pairing interaction 
are likely to reduce d-wave superconductivity. 

Since we find T c oc J e f / , it may be possible to increase 
T c by enhancing J e f f . In our model we find that a weak 
inhomogeneity decreases Jeff, but other kinds of inhomo- 
geneity might increase J e / / and thus, if our assumption is 
correct, enhance T c . For example, weak on-site disorder 
seems to increase J in the vicinity of defects [HI . 

Regarding the speculation for the existence of an opti- 
mal inhomogeneity raised in Ref. [l7j |. our results based 
on the N c = 4 cluster indicate that T c is an increasing 
function of t' in the weakly inhomogeneous system and 
the optimal t' occurs at t' = t. Therefore our results for 
the weakly inhomogeneous system are not in agreement 
with this speculation , but they do not exclude the possi- 
bility for an optimal t' at smaller values of t' . In addition, 
it will also be interesting to see how T c scales with cluster 
size N c when the weak inhomogeneity is present. 

Conclusion. Using the Dynamical Cluster Approxima- 
tion for an N c = 4 cluster we find that a weak inhomo- 
geneity in the checkerboard Hubbard model suppresses 
d-wave superconductivity. The characteristic spin exci- 
tation energy and the strength of the pairing interaction 
decrease along with decreasing T c suggesting a strong 
correlation between these quantities. 
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